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ABSTRACT 


This  study  examines  earthquakes  in  Kamchatka,  the  Tien  Shan,  the  Pamirs, 
the  Baikal  rift  zone,  and  the  Caucasus  and  explosions  in  the  Aleutians,  East 
Kazakh,  West  Kazakh,  Southwest  Russia,  Lop  Nor,  and  Baikal  in  a  seismic  dis¬ 
crimination  context.  Because  of  their  low  magnitudes  (m^  from  4.2  to  6.0 

and  M  from  3.2  to  5.8)  and  their  varying  propagation  effects,  it  was  impos- 

s 

sible  on  the  basis  of  first  motions  or  LR  amplitudes  to  determine  source 
mechanisms  for  most  of  the  earthquakes.  Seismograms  from  the  arrays  ALPA, 
LASA,  and  NORSAR,  the  HGLP-  and  the  WWSSN  stations  were  analyzed  for  m^,  Mg, 
corner  frequency,  long-period  spectral  level,  long-period  and  short-period 
body-wave  excitation,  complexities,  spectral  ratios,  radiation  patterns, 
depth  of  focus,  and  higher-mode  surface  waves.  The  major  successful  discrim¬ 
inants  for  our  data  set  were  found  to  be  M^-m^  and  identification  of  pP.  In 
a  stepwise  linear  discrimination  analysis  using  the  discrimination  parameters 
n^,  Ms>  corner  frequency,  complexity,  spectral  ratio,  and  long-period  spec¬ 
tral  level,  the  earthquakes  were  separated  from  the  explosions  with  high  con¬ 
fidence.  In  addition,  some  separation  by  region  was  evident  for  the  earth¬ 
quakes  alone. 

However,  it  is  clear  that  in  many  regions  a  shot  array  could  be  used  to 
successfully  evade  discrimination  which  would  then  have  to  be  made  using  other 
means  such  as  depth,  relative  P  and  pP  amplitudes,  and  ratio  of  compressional 
to  shear  phases  at  similar  frequencies. 
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INTRODUCTION 


Variation  of  source  characteristics  and  path  effects  from  events  in  vari¬ 
ous  regions  is  likely  to  compound  the  problem  of  seismic  discrimination  of 
explosions  and  earthquakes,  unless  reasonably  accurate  descriptions  of  the 
cause  and  degree  of  such  variations  are  obtained.  To  accomplish  this  the 
authors  have  recently  attempted,  for  five  different  source  regions  in  and 
around  the  Union  of  Soviet  Socialist  Republic  (USSR)  (von  Seggern  and  Sobel, 
1976;  Sobel  and  von  Seggern,  1976;  Sobel  et  al.,  1977a,  1977b,  and  1977c),  to 
identify  the  major  source  and  path  effects  which  impact  upon  discrimination  of 
Asian  events  and  to  assess  likelihood  of  correct  event  identification  in  these 
regions.  However,  the  problem  of  regional  effects  can  be  circumvented  by  lim¬ 
iting  the  computation  of  a  discriminant  function  for  earthquakes  and  explosions 
to  events  within  a  geophysically  distinct  region.  This  type  of  analysis  was 
conducted  on  a  suite  of  events  in  the  Southwestern  United  States  (von  Seggern 
and  Rivers,  1977)  with  encouraging  results  for  discrimination  on  a  regional 
basis. 


von  Seggern,  D.  H.,  and  P.  A.  Sobel,  1976.  Study  of  selected  Kamchatka  earth¬ 
quakes  in  a  seismic  discrimination  context,  SDAC  Report  No.  TR-76-10, 
Teledyne  Geotech,  Alexandria,  Va.  22314. 

Sobel,  P.  A.,  and  D.  H.  von  Seggern,  1976.  Study  of  selected  events  in  the 
Tien  Shan  Region  in  a  seismic  discrimination  context,  SDAC  Report  No. 
TR-76-9,  Teledyne  Geotech  Alexandria,  Va.  22314. 

Sobel,  P.  A.,  D.  H.  von  Seggern,  E.  I.  Sweetser,  and  D.  W.  Rivers,  1977a. 

Study  of  selected  events  in  the  Pamirs  in  a  seismic  discrimination  con¬ 
text,  SDAC  Report  No.  TR-77-3,  Teledyne  Geotech,  Alexandria,  Va.  22314. 

Sobel,  P.  A.,  D.  H.  von  Seggern,  E.  1.  Sweetser,  and  D.  W.  Rivers,  1977b. 

Study  of  selected  events  in  the  Baikal  rift  zone  in  a  seismic  discrimin¬ 

ation  context,  SDAC  Report  No.  TR-77-5,  Teledyne  Geotech,  Alexandria, 

Va.  22314. 

Sobel,  P.  A.,  D.  H.  von  Seggern,  E.  I.  Sweetser,  and  D.  W.  Rivers,  1977c. 

Study  of  selected  events  in  the  Caucasus  in  a  seismic  discrimination  con¬ 

text,  SDAC  Report  No.  TR-77-6,  Teledyne  Geotech,  Alexandria,  Va.  22314. 

von  Seggern  D.  H. ,  and  W.  D.  Rivers,  1977.  Seismic  discrimination  of  earth¬ 
quakes  and  explosions  with  applications  to  the  Southwestern  United 
States,  SDAC  Report  No.  TR-77-10,  Teledyne  Geotech,  Alexandria,  Va. 

22314. 
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The  purpose  of  this  report  is  to  synthesize  and  analyze  the  work  contained 
in  the  five  regional  Asian  studies.  Through  this  presentation,  the  diversi¬ 
ties  and  similarities  of  the  recorded  data  can  be  appreciated  because  a 
nearly  uniform  approach  to  the  study  of  each  of  the  five  regions  was  made. 

In  this  report,  the  five  previous  efforts  are  collectively  labeled  the 
"Asian  regional  studies"  or  "Asian  regional  reports." 

The  scheme  of  this  report  is  as  follows:  First,  because  the  basic  illus¬ 
trations  of  the  signal  waveforms  are  one  of  the  more  important  parts  of  this 
document,  the  data  used  in  the  Asian  regional  studies  will  be  reviewed  and 
depicted.  For  short-period  data,  primarily  LASA  and  NORSAR  recordings  were 
used,  with  some  supplementary  work  on  WWSSN  recordings  for  first  motions  and 
magnitude.  For  long -period  data,  the  large  arrays  ALPA,  LASA,  and  NORSAR  were 
used,  with  additional  readings  from  the  HGLP  network  and  the  WWSSN  network  for 
magnitude.  The  dates  of  events  in  the  five  Asian  regional  studies  were  from 
1971  to  1975,  before  the  advent  of  the  SRO  network  of  high-quality  stations. 

Second,  the  signal  analysis  sections  of  the  Asian  regional  reports  are  sum¬ 
marized,  leading  to  discussions  of  various  discrimination  parameters  and  the 
effects  of  source  mechanisms  and  travel  path  on  them.  Discrimination  para¬ 
meters  for  all  five  regions  will  be  plotted  together  without  "regional"  cor¬ 
rections.  The  discriminants  treated  are:  Ms,  m^,  corner  frequency,  long- 
period  spectral  level,  complexity,  spectral  ratio,  radiation  pattern,  depth  of 
focus,  higher-mode  generation,  and  shear-waves. 

Finally,  new  work  is  added  on  stepwise  classification  experiments  not  con- 
tained  in  the  original  Asian  regional  reports.  This  section  will  treat  clas¬ 
sification  of  events  on  a  regional  basis  and  classification  of  all  the  Asian 
events  together.  In  this  manner,  the  relative  value  of  various  discriminant 
parameters  for  these  events  will  be  established. 
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TECTONIC  SETTING 


The  earthquakes  studied  here  occurred  in  five  different  tectonic  settings: 
the  Kamchatka  Peninsula,  the  Tien  Shan  uplift,  the  Pamirs  uplift,  the  Baikal 
rift  zone,  and  the  Caucasus  fold  belt.  The  locations  of  these  earthquakes  are 
depicted  on  Figure  1.  Each  of  the  regions  is  seismic;  and  they  have  therefore 
been  extensively  studied  by  Soviet  seismologists  (Medvedev,  1968).  The  Tien 
Shan,  Pamir,  and  Caucasus  earthquakes  lie  in  distinct  seismic  regions,  but  all 
are  undergoing  uplift  as  a  result  of  the  collision  of  the  Indian  and  Eurasian 
plates  (Molnar  and  Tapponier,  1975).  The  Baikal  rift  zone  is  a  region  of 
crustal  extension  which  Molnar  and  Tapponier  thought  to  be  a  distant  manifes¬ 
tation  of  the  stress  system  created  by  the  Indian-Eurasian  collision.  The 
Kamchatka  seismic  region  is  unique  in  this  study  because  it  is  connected  with 
the  underthrusting  of  the  Pacific  plate. 

The  Kamchatka  earthquakes  occurred  along  the  coastline  within  the  upper 
margin  of  the  dipping  Pacific  plate  (25-50  km  deep).  A  typical  source  mechan¬ 
ism  for  these  events  would  be  a  fault  plane  striking  NNE  parallel  to  the  tec¬ 
tonic  trend  and  dipping  approximately  25  degrees  WNW.  The  events  are  charac¬ 
teristically  thrust  faults  with  a  slip  angle  of  roughly  90  degrees.  The  maxima 
of  the  P  radiation  pattern  is  downward,  giving  positive  first  motions  to  all 
teleseismic  signals.  Most  of  the  earthquakes  in  the  Tien  Shan  area  are  asso¬ 
ciated  with  large,  sub-parallel,  strike-slip  faults  inferred  from  analysis  of  a 
satellite  photomosaic  of  the  area.  Review  of  earthquake  focal  mechanisms 
in  the  Tien  Shan  show  both  strike-slip  and  thrust  movements  mainly  along  planes 
with  steep  dip  angles  (Sobel  et  al.,  1977a).  The  seismicity  in  the  Caucasus 
Mountains  generally  follows  the  northwest-southeast  trend  of  the  principal 
geological  features:  thrust  faulting,  with  the  fault  planes  perpendicular  to 
the  trend  of  the  geological  structures,  is  dominant,  but  considerable  variance 
to  the  fault-plane  orientations  exists  and  strike-slip  motion  is  not  uncommon 
(Sobel  et  al.,  1977c).  In  the  Baikal  rift  zone,  normal  faulting  is  prevalent; 


Molnar,  P.,  and  P.  Tapponnier,  1975.  Cenozoic  tectonics  of  Asia:  effects  of 
a  continental  collision.  Science.  189,  419. 

Medvedev,  S.  V.  (ed.),  1968.  Seismic  Zoning  in  the  USSR,  Keter  Publ.  House, 
Jerusalem,  Israel  (English  translation). 
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EARTHQUAKES 


Map  of  Eurasia  showing  the  earthquakes  and  explosions  used  in 
this  study. 


and  faults  bordering  a  system  of  grabens  are  clearly  observed  on  satellite 
photographs  of  the  area  (Sobel  et  al.,  1977b).  With  normal  faulting  the  nodal 
plane  of  the  P  radiation  pattern  is  downward,  giving  positive  and  negative  first 
motions  to  teleseismic  signals. 
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DATA 


Earthquake  Selection 

Since  the  emphasis  is  on  discrimination,  only  crustal  depth  earthquakes, 
except  one  subcrustal  event  in  the  Caspian  Sea,  were  chosen  for  this  study. 
Earthquakes  were  selected  in  the  years  1971-1975  so  digital  data  could  be  uti¬ 
lized  from  the  ALPA,  LASA,  and  NORSAR  arrays  and  the  HGLP  network.  In  the 
Kamchatka,  Tien  Shan,  and  Pamir  regions,  earthquakes  were  chosen  that  had  low 
reported  Mg  for  their  so  that  these  events  fall  close  to  the  explosion  pop¬ 
ulation  on  an  M^m^  diagram.  In  the  Baikal  and  Caucasus  regions  such  low  Mg 
events  were  not  found,  and  all  crustal  events  were  selected  in  the  years  1971- 
1975  for  which  array  data  could  be  recovered  with  at  least  fair  signal-to-noise 
ratios.  These  few  earthquakes  cannot  be  claimed  to  represent  adequately  the 
true  variability  within  the  five  source  regions,  but  they  are  representative 
of  the  type  of  earthquakes  most  troublesome  in  a  discrimination  context.  The 
epicentral  data  for  the  earthquakes  is  listed  in  Table  I,  and  the  locations 
are  shown  in  Figure  1. 

Explosion  Selection 

In  the  five  Eurasian  reports,  earthquakes  in  the  regions  of  study  were 
compared  to  explosions  detonated  as  closely  as  possible  to  these  seismic 
regions.  Explosions  were  selected  from  the  Aleutians  (MILROW),  East  Kazakh, 
West  Kazakh,  Southwest  Russia,  Lop  Nor,  and  Baikal.  The  explosion  epicenter 
data  and  the  earthquake  regions  to  which  they  are  compared  are  listed  in 
Table  II,  and  their  locations  are  shown  in  Figure  1. 

Seismic  Stations 

Digital  data  from  the  arrays  ALPA,  LASA,  and  NORSAR  and  from  the  available 
HGLP  stations  were  gathered  for  these  events.  A  significant  percentage  of  the 
HGLP  data  either  could  not  be  recovered  or  it  was  not  usable.  Also  examined 
were  short-period  and  long-period  data  on  film  chips  from  reliable  sites  in 
the  WWSSN  to  supplement  the  magnitude  calculations  and  focal-mechanism  investi¬ 
gations.  WWSSN  stations  were  selected  on  the  basis  of  magnification,  proximity 
to  the  regions  of  study,  and  azimuthal  relation  to  the  source  regions.  The 
Baikal  explosion  was  added  late  to  the  data  set;  and  since  WWSSN  data  of  the 
Baikal  explosion  was  not  available  at  the  time,  SRO  data  for  this  event  was 
used  instead. 
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Bungum  and  Tjostheim,  1976 


Short-period  Waveforms  at  LAS A  and  NORSAR 

Short-period  P  phases  recorded  at  LASA  and  NORSAR  for  the  five  earthquake 
regions  are  shown  in  Figures  2  through  6;  those  for  the  explosions  are  shown 
in  Figure  7.  Phased  AO  subarray  beams  from  LASA  and  phased  C3  subarray  beams 
from  NORSAR  were  used  to  avoid  signal  loss  at  high  frequencies  due  to  full- 
array  beaming.  Where  there  are  no  records  shown,  the  data  for  that  time 
period  was  unrecoverable  or  there  was  no  detected  signal. 

The  short-period  P  waveforms  for  the  Kamchatka  earthquakes  are  shown  in 
Figure  2.  Although  a  source  mechanism  solution  exists  for  only  one  of  these 
events,  the  730304  (year-month-day)  earthquake,  the  association  of  these  events 
with  the  under thrusting  of  the  oceanic  plate  and  the  dominance  of  teleseismic 
compressional  first  motion  indicate  that  these  are  thrust  faults.  In  addition, 
synthetic  seismograms  computed  for  the  730304  earthquake  match  the  directions 
of  first  motion,  arrival  times  of  pP,  P-wave  amplitudes,  and  the  shape  of  the 
first  few  cycles  of  observed  long-period  P  at  LASA,  NORSAR,  and  eight  world¬ 
wide  stations  (von  Seggern  and  Sobel,  1976).  The  pP  phase  is  observed  on  all 
LASA  records  but  only  on  two  of  nine  NORSAR  records.  Observations  of  pP  at 
NORSAR  are  lacking  because  the  pP  raypath  to  the  surface  for  NORSAR  lies  near 
the  auxiliary  focal  plane,  where  the  amplitudes  are  low.  This  difference  in 
the  P  to  pP  ratios  serves  as  a  powerful  counter  evasion  technique  to  shot 
arrays.  For  LASA,  the  predicted  amplitude  of  pP  is  close  to  that  of  P,  and 
this  is  observed  in  most  cases.  There  were  similarities  among  many  of  the 
initial  P  waveforms  at  LASA,  but  the  waveforms  recorded  at  NORSAR  show  little 
uniformity  and  are  always  of  higher  frequency  content  than  corresponding  LASA 
signals. 

Short-period  P  waveforms  for  the  Tien  Shan  earthquakes  are  shown  in  Figure 

3.  Definite  source  mechanisms  could  not  be  determined  for  any  of  these  events 

using  teleseismic  data  due  to  their  small  size  (m,  <  5.0);  satellite  photo- 

b  — 

graphs  suggest  that  the  mechanisms  could  be  either  thrust  faulting  or  strike- 
slip  faulting  (Sobel  and  von  Seggern,  1976).  There  are  many  differences  in 
signal  shape  and  frequency  content  of  these  P  waveforms  at  both  LASA  and 
NORSAR,  and  again  the  NORSAR  signals  have  higher  frequency  content  than  the 
LASA  ones.  Clearly,  no  one  waveform  could  be  used  to  approximate  the  P  sig¬ 
nals  at  either  LASA  or  NORSAR.  The  pP  phases  were  not  well  recorded  at  either 
LASA  or  NORSAR  for  the  Tien  Shan  events. 
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Figure  3  LASA  and  NORSAR  short-period  P  waveforms  for  the  ' 
earthquakes. 
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LASA  and  NORSAR  short-period  waveforms  for  the  Pamir  earthquakes 
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Figure  6  LASA  and  NORSAR  short-period  P  waveforms  for  the  Caucasus 
earthquakes . 
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Figure  7  (Cont.)  LASA  and  NORSAK  short-period  P  waveforms  for  the 
explosions  used  in  this  study. 


Short-period  P  waveforms  for  the  Pamir  earthquakes  are  shown  in  Figure  4. 
Again,  focal  mechanism  solutions  for  these  events  could  not  be  determined  due 
to  their  small  size  (m^  <^4.8).  For  these  events  the  presence  of  negative  and 
positive  first  motions  at  teleseismic  distances,  without  any  definite  overall 
pattern,  indicates  that  they  are  not  thrust  faults  and  that  there  is  probably 
no  common  source  mechanism  for  them  (Sobel  et  al.,  1977a).  The  phase  pP  was 
present  on  some  of  the  LASA  and  NORSAR  records.  There  are  many  differences  in 
the  sign  of  first  motion,  frequency  content,  and  waveform  among  these  events; 
such  dissimilarity  is  perhaps  unexpected  because  of  the  epicenters  of  all  but 
the  730103  and  the  740220  events  were  located  within  a  one-degree  square  on 
the  earth's  surface.  Note,  however,  that  this  degree  of  variability  was  also 
found  among  recordings  of  Kamchatka  events  within  such  a  small  area  by  Bland- 
ford  and  Clark  (1975).  Clearly,  no  one  signal  represents  the  waveforms  at 
LASA  or  NORSAR  for  these  Pamir  events.  The  frequency  content  of  the  NORSAR 
signals  is  strikingly  higher  than  that  of  the  LASA  ones. 

Short-period  P  waveforms  for  the  Baikal  earthquakes  are  shown  in  Figure 
5.  Again,  small  magnitudes  (m^  <_  5.3)  preclude  determination  of  source  mechan¬ 
isms  for  these  events.  However,  the  predominance  of  negative  first  motions  at 
teleseismic  distances  suggests  that  these  events  are  normal  faults  on  steeply 
dipping  planes,  which  agrees  with  known  focal  mechanisms  for  this  region  (Sobel 
et  al.,  1977b).  Note  that  some  LASA  and  NORSAR  first  motions  are  positive, 
indicating  the  vagaries  of  presumed  first-motion  directions  for  small-magni¬ 
tude  events.  None  of  the  waveforms  are  similar,  and  the  contrast  between 
NORSAR  and  LASA  recordings  is  pronounced  for  at  least  two  of  the  events, 

710614  and  740621.  The  pP  was  observed  at  LASA  or  NORSAR  for  only  three  of 
the  events,  but  in  each  case  identification  is  not  certain  and  resulted  only 
from  corroborating  WWSSN  recordings. 

Short-period  LASA  and  NORSAR  P  waveforms  for  the  Caucasus  earthquakes  are 
shown  in  Figure  6.  Again,  we  were  not  able  to  determine  the  focal  mechanisms 
due  to  the  small  size  (m^  <_  5.1)  of  the  events.  The  variable  patterns  of 

Blandford,  R. ,  and  D.  Clark,  1975.  Variability  of  seismic  waveforms  recorded 

at  LASA  from  small  subregions  of  Kamchatka,  SDAC  TR-75-12,  Teledyne  Geo¬ 
tech,  Alexandria,  Va.  22314. 


negative  and  positive  first  motions  at  teleseismic  distances  suggest  that  these 
are  not  thrust  faults  and  that  probably  no  common  source  mechanism  for  all 
these  events  exists  (Sobel  et  al.,  1977c).  Waveforms  for  the  events  were  not 
similar  at  either  LASA  or  NORSAR.  The  phase  pP  was  apparently  recorded,  al¬ 
though  poorly  in  most  cases,  at  either  LASA  or  NORSAR  for  all  the  earthquakes. 
In  contrast  to  the  previous  four  regions,  the  NORSAR  signals  for  these  Caucas¬ 
us  events  do  not  appear  to  have  higher  frequency  content  than  LASA  recordings. 

Short-period  P  waveforms  for  the  explosions  used  in  this  study  are  shown 
in  Figure  7.  For  three  of  the  explosions,  no  signals  were  recorded  at  LASA 
due  to  the  small  size  of  the  events  (m^  _<  4.8).  Although  several  waveforms 
from  nearby  events  were  present  at  LASA  and  NORSAR  from  the  two  Kazakh  test 
sites,  only  the  East  Kazakh  signals  recorded  at  LASA  showed  enough  similarity 
that  one  waveform  represented  the  entire  group.  This  is  true  in  spite  of 
the  fact  that  the  explosion  locations  within  either  East  or  West  Kazakh 
(Table  II)  were  in  close  proxim. 

Long-period  waveforms  at  NORSAR 

Of  the  three  long-period  arrays  used,  only  NORSAR  had  enough  good  surface 
waves  recorded  to  warrant  discussion  here.  All  observed  Rayleigh  waves  at 
NORSAR  for  the  earthquake  regions  are  shown  in  Figures  8  through  12;  those 
for  the  explosions  are  shown  in  Figure  13.  Missing  records  for  an  event 
indicate  that  no  data  or  no  detected  signal  for  the  time  period  existed. 

The  beams  are  all  formed  with  delays  corresponding  to  a  velocity  of  3.6  km/sec. 
Recordings  displayed  here  begin  at  approximately  the  group  arrival  time  of  4.5 
km/sec  energy  and  end  at  the  group  arrival  time  of  2.7  km/sec  energy,  and  they 
are  roughly  aligned  according  to  time  after  origin  time. 

Figure  8  shows  the  long-period  NORSAR  waveforms  for  the  Kamchatka  earth¬ 
quakes.  Epicentral  distance  is  roughly  65°.  All  the  traces  show  effects  of 
lateral  refractions  predicted  by  ray  tracing  for  20-second  Rayleigh  waves  from 
Kamchatka  to  NORSAR  in  von  Seggern  and  Sobel  (1976).  The  maxima  in  the  period 
range  17-23  sec,  normally  used  for  Mg  estimation,  arrives  with  a  group  velocity 
of  roughly  3.0  km/sec  in  all  cases.  While  there  was  some  similarity  among  the 
waveforms,  there  was  also  considerable  variability  relative  to  the  small 
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NORSAR  long-period  waveforms  for  the  Kamchatka  earthquakes 
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NORSAR  long-period  waveforms  for  the  Tien  Shan  earthquakes 
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Figure  10  NORSAR  long-period  waveforms  for  the  Pamir  earthquakes 
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Figure  11  NORSAR  long-period  waveforms  for  the  Baikal  earthquakes. 
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Figure  12  NORSAR  long-period  waveforms  for  the  Caucasus  earthquakes. 
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NORSAR  long-period  waveforms  for  representative  explosions 


dimension  of  the  source  area.  The  least  similar  record  of  the  group  is  for  the 
740427  event,  which  is  also  the  farthest  north  of  the  five  events.  The  ray 
paths  from  Kamchatka  to  NORSAR  traverse  the  continental  margins  of  Northeast 
Asia  and  Northern  Scandinavia,  both  sites  of  possible  lateral  refractions  which 
could  create  the  multipathing  observed  at  NORSAR.  The  waveforms  may  also  be 
affected  by  the  spatial  relation  of  the  exact  hypocenters  to  the  structure  of 
the  dipping  lithospheric  plate  under  Kamchatka  and  by  the  frequency-dependent 
source  excitation. 

Figure  9  shows  the  NORSAR  long-period  waveforms  for  the  Tien  Shan  earth¬ 
quakes.  The  epicentral  distance  is  roughly  45°.  All  records  have  approxi¬ 
mately  the  same  arrival  time  of  20-second  energy,  contained  within  a  distinct 
Airy  phase.  Ray  paths  to  NORSAR  are  entirely  continental,  which  explains  the 
arrival  time  and  shape  of  this  Airy  phase.  No  evidence  exists  of  multipathing 
predicted  by  the  20-second  Rayleigh-wave  ray  tracing  of  Sobel  and  von  Seggern 
(1976),  but  the  predicted  phase  arrival  times  were  not  greatly  different. 

While  the  first  five  events  have  several  cycles  in  the  Airy  phase,  and  the 
waveforms  show  some  similarity,  events  750114  and  750222  apparently  Lave  a 
more  concentrated  Airy  phase.  These  latter  two  events  occurred  geographi¬ 
cally  in  the  middle  of  the  other  five  events. 

The  NORSAR  long-period  records  for  the  Pamir  events  are  shown  in  Figure 
10.  The  epicentral  distance  is  roughly  40°.  Little,  if  any,  multipathing  is 
present,  and  the  low  Rayleigh-wave  group  velocity  for  20-sec  energy  reflects 
the  continental  path  of  the  rays  from  the  Pamirs  to  NORSAR.  An  Airy  phase  is 
distinct  only  on  some  recordings,  and  it  is  not  as  strong  as  for  the  Tien  Shan 
events.  All  except  the  first  four  events  occur  within  a  small  geographic 
area,  and  yet  the  waveforms  recorded  at  NORSAR  show  only  some  similarity. 
Although  differences  between  the  second  and  third  740811  events,  for  instance, 
could  be  attributed  to  focal  depth,  our  depth  determinations  (Table  I)  indi¬ 
cated  that  they  had  similar  depths.  Tectonically  active  structures  in  the 
Northern  Pamirs  may  in  part  be  responsible  for  differences  in  waveforms  from 
event  located  in  this  region. 

Figure  11  shows  the  NORSAR  long-period  signals  for  events  in  the  Baikal 
rift  zone.  The  epicentral  distances  are  roughly  50°.  There  is  evidence  of 
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multipathing  on  the  720115  and  721125  records.  All  the  waveforms  are  differ¬ 
ent  because  of  refraction  effects  within  the  Baikal  rift  zone  or  because  the 
rays  pass  through  the  Ural  Mountains.  The  720809  waveform  is  especially  dif¬ 
ferent,  and  the  location  of  this  event  (Table  I)  is  apart  from  others  of  the 
group.  The  late  arrivals  of  20-second  energy  (2.9  to  3.1  km/sec)  reflect  the 
totally  continental  ray  paths. 

NORSAR  long-period  records  for  the  Caucasus  events  are  shown  in  Figure 
12.  The  epicentral  distance  is  roughly  30°.  The  great  dissimilarity  among 
these  waveforms  may  partially  result  from  refraction  effects  as  they  pass 
through  the  northern  border  of  the  Black  Sea  and  Caspian  Sea  sedimentary 
depression.  Note  that  while  events  711015  and  731214  are  located  close  geo¬ 
graphically,  they  have  totally  dissimilar  waveforms.  This  lack  of  shorter- 
period  energy  for  the  731214  event  is  supporting  evidence  for  the  reported 
NEIS  subcrustal  depth  of  79  km.  The  Caucasus  waveforms  do  not  show  distinct 
Airy  phases  such  as  those  seen  on  the  Tien  Shan  or  Pamir  waveforms. 

Figure  13  shows  the  Rayleigh-wave  signals  recorded  at  NORSAR  for  two 
East  Kazakh  explosions  (A  'v  38°),  located  in  close  geographical  proximity  and 
two  Southwest  Russia  explosions  (A  'v  25°),  separated  by  approximately  three 
degrees  of  epicentral  distance.  For  each  pair  the  waveforms  are  similar.  The 
low  Rayleigh-wave  velocities  reflect  the  continental  ray  paths  to  NORSAR.  No 
multipathing  is  observed  in  these  records.  Note  that  NORSAR  did  not  detect 
Rayleigh  waves  from  either  the  Baikal  or  the  Lop  Nor  explosion  in  Table  II. 
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SIGNAL  ANALYSIS 


In  this  section  are  discussed  effects  of  source  radiation  pattern,  source 
spectral  shape,  and  inelastic  attenuation  on  observed  body-wave  and  surface- 
wave  signals  for  events  in  this  study. 

Source  Mechanisms 

(  Efforts  were  made  to  identify  the  effects  of  source  mechanism  on  observed 

body-wave  and  surface-wave  phases  from  the  earthquakes  studied  here.  From  a 
knowledge  of  the  tectonic  settings  of  these  earthquake,  the  Kamchatka  events 
are  assumed  to  be  thrust  faults,  the  Tien  Shan,  Pamir,  and  Caucasus  events 
are  assumed  to  be  strike-slip  or  thrust  faults,  and  the  Baikal  events  are 
assumed  to  be  normal  faults.  The  only  well-defined  source  mechanism  for  any 
of  the  earthquakes  studied  here  was  for  the  730304  event  in  Kamchatka  (Stauder 
and  Maulchin,  1976),  a  thrust  fault  and  the  largest  event  in  our  earthquake 
data  set  (m^  =  6.0).  Determination  of  the  fault  plane  for  an  earthquake 
whose  magnitude  is  below  6.0  is  generally  not  reliable  with  teleseismic  data 
because  P  waves  are  not  recorded  on  the  long-period  seismograms  and,  except 
for  the  730304  Kamchatka  event,  all  earthquakes  have  <  6.  Short-period 
first  motions  have  been  read  from  array,  HGLP ,  and  WWSSN  recordings  for  all 
!  the  earthquakes  in  the  study.  An  example  typical  of  the  ambiguous  data  sets 

i  for  our  events  is  shown  in  Figure  14.  On  the  basis  of  such  plots  alone, 

I  fault  planes  could  not  be  determined  because  dilatational  and  compressional 

J  first  motions  did  not  separate.  A  composite  plot  of  all  Kamchatka  data  showed 

a  predominance  of  positive  first  motions,  suggesting  that  these  events  are 
indeed  thrust  faults.  A  similar  composite  plot  of  all  the  Baikal  data  showed 
a  predominance  of  negative  first  motions,  suggesting  that  the  Baikal  events 
are  normal  faults.  No  common  source  mechanism  was  clearly  indicated  for  the 
j  Tien  Shan,  Pamir,  or  Causasus  earthquakes. 

i  Attempts  were  made  to  find  the  LR  amplitude  radiation  pattern  for  sev- 

1 , 

eral  of  the  events  by  plotting  the  antilogs  of  the  Mg  values.  Typical  results 
for  an  event  are  shown  in  Figure  15.  (The  maximum  antilog  is  plotted  on  the 

i.(  Stauder,  W. ,  and  L.  Mualchin,  1976.  Fault  motion  in  the  large  earthquakes  of 

ij  the  Kurile-Kamchatka  Arc  and  of  the  Kurile-Hokkaido  corner,  J.  Geophys. 

!'  Res. ,  81,  297. 
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outer  circle,  and  all  other  antilog  values  are  scaled  linearly  between  the 
center  of  the  plot  and  the  radius  of  the  circle.  (The  data  point  with  an  arrow 
is  a  noise  measurement).  Few  observed  distributions  of  amplitudes  conformed 
to  any  quadripole  radiation  pattern.  The  scatter  of  data  points  was  a  result 
of  the  large  epicentral  distances  involved,  because  deviations  due  to  propa¬ 
gation  path  differences  accumulated  to  the  point  where  they  overwhelmed  the 
initial  radiation  pattern  effect  (Aki,  1970).  For  Rayleigh  waves  of  20-sec¬ 
onds  period,  a  method  of  predicting  amplitudes  by  geometrical  ray  tracing  has 
been  presented  by  von  Seggern  and  Sobel  (1975).  Figure  16  illustrates  the 
path  of  20-second  Rayleigh  waves  that  emanated  from  an  event  in  the  Tien  Shan 
Mountains.  The  refraction  effect  of  the  earth's  non-uniformity  in  20-second 
LR  phase  velocity  causes  focusing,  defocusing,  and  multipathing  and,  therefore, 
large  teleseismic  amplitude  variations.  Mg  residuals  can  be  predicted  in  some 
cases  from  the  ray-pattern  density  near  each  of  the  stations.  No  attempt  was 
made  to  relate  ray-tracing  predictions  to  observed  LR  amplitudes  since  we 
could  not  correct  for  source  mechanism  and  did  not  have  enough  information  to 
correct  for  lateral  differences  in  attenuation.  Source  locations  only  100  km 
apart  on  the  earth's  surface  will  produce  large  amplitude  differences  at  cer¬ 
tain  stations  due  to  changing  patterns  of  focusing  and  defocusing  (von  Seggern 
and  Sobel,  1976).  This  phenomenon  precludes  determining  source  radiation  pat¬ 
terns  using  a  master-event  technique,  unless  there  is  a  small  separation 
between  epicenters. 

Corner  Frequencies,  Long-Period  Spectral  Levels,  and  Seismic  Moments 

In  the  five  Eurasian  reports,  amplitude  spectra  were  displayed  for  the 
LASA  and  NORSAR  short-period  P  recordings  for  all  the  earthquakes  and  explo¬ 
sions.  Some  examples  of  these  spectra  are  shown  in  Figure  17  to  illustrate 
how  spectral  parameters  were  determined.  The  mean  of  the  noise  sample  was 
removed  from  the  signal  and  noise  samples  shown  on  the  left  side  of  each 
figure.  Then  the  signal  and  noise  samples  were  tapered  by  one  half  a  cosine 
bell  over  6  pts  at  each  end  of  the  6.4  second  windows  (10  s/sec).  Instrument 


von  Seggern,  D.  H. ,  and  P.  Sobel,  1975.  Experiments  in  refining  M  estimates 
for  seismic  events,  SDAC  Report  No.  TR-75-17,  Teledyne  Geotecli,  Alexan¬ 
dria,  Virginia  22314. 
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Figure  17  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P  waves  from  the 

740226  Kamchatka  earthquake,  the  720102  Tien  Shan  earthquake 
and  the  730815  W.  Kazakh  explosion. 
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Figure  17  (cont.)  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P  waves  from 

the  740226  Kamchatka  earthquake,  and  720102  Tien  Shan 
earthquake,  and  the  730815  W.  Kazakh  explosion. 
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Figure  17  (cont.)  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P  waves  from 

the  740226  Kamchatka  earthquake,  and  720102  Tien  Shan 
earthquake,  and  the  730815  W.  Kazakh  explosion. 
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LASA  AO  and  NORSAR  C3  subarray  spectra  of  P  waves  from 
the  740226  Kamchatka  earthquake,  and  720102  Tien  Shan 
earthquake,  and  the  730815  W.  Kazakh  explosion. 


Figure  17  (cont.)  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P  waves  from 

the  740226  Kamchatka  earthquake,  and  720102  Tien  Shan 
earthquake,  and  the  730815  W.  Kazakh  explosion. 


response  was  removed  from  the  raw  spectra  to  get  ground  displacement  (t*=0), 
but  noise  spectra  have  not  been  subtracted  from  the  signal  spectra.  Attenua¬ 
tion  was  removed  by  multiplying  by  the  factor  exp  [nft*]  with  t*  values  (deter¬ 
mined  in  the  next  section  of  this  report)  in  order  to  correct  them  back  to  the 

-2 

source.  An  f  source  model  was  used  for  both  earthquakes  and  explosions, 

-3 

and  an  f  source  model  was  also  used  for  earthquakes  only. 

Corner  frequencies,  long-period  spectral  levels,  and  seismic  moments 

have  been  estimated  from  the  resulting  source  spectra.  Corner  frequencies 

-2 

were  estimated  with  the  assumptions  of  complete  stress  drop  and  an  f  or 
-3 

f  asymptotic  relation  at  high  frequencies.  Both  corner  frequencies  (f  )  and 

A 

long-period  spectral  levels  (Uq)  are  shown  in  Figure  17.  Occasionally  spectra 
for  25.6-second  windows  were  computed  for  cases  where  the  signal-to-noise 
ratio  was  high.  However,  for  these  cases  the  shapes  of  the  spectra  were  suf¬ 
ficiently  close  to  those  for  the  6.4-second  windows  that  little  or  no  addi¬ 
tional  information  was  gained. 


Seismic  moment  was  calculated  from  the  long-period  spectral  levels 
| U0 1  using  the  relation  (Hanks  and  Thatcher,  1972) 


M  = 
o 


4irpa3  ju 


%D 


where  R_,  was  assumed  to  be  unity  due  to  lack  of  knowledge  of  the  focal  mech- 
t)<p 

anism,  values  of  density  (p)  and  compressional  velocity  (a)  appropriate  to 
the  source  region  and  depth  were  used,  and  the  divergence  factor  D  was  applied 
(Ben-Menahem  et  al.,  1965).  Graphs  of  moment  versus  corner  frequency  are 
shown  in  Figures  18  and  19  for  all  earthquake  data  at  LASA  and  NORSAR.  The 
events  fall  between  or  close  to  the  100-bar  and  1-bar  stress  drop  lines  of 
Hanks  and  Thatcher  (1972),  indicating  that  they  are  generally  of  intermediate 
stress  drop.  The  Kamchatka  events,  which  are  generally  of  higher  magnitudes 
than  earthquakes  from  the  other  regions,  generally  possess  higher  moments  and 
greater  stress  drops.  The  Baikal  events  have  generally  lower  stress  drops 
than  the  other  earthquakes.  The  Mq  vs.  f  plots  can  be  interpreted  in  terms 


Hanks,  T. ,  and  W.  Thatcher,  1972.  A  graphical  representation  of  seismic 
source  parameters,  J.  Geophy.  Res.,  77 ,  4393. 
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moment  versus  corner  frequency  for  earthquake 
from  LASA  P  recordings. 


NORSflR  DfiTR 


Seismic  moment  versus  corner  frequency  for  earthquakes  from  al] 
regions  from  NORSAR  P  recordings. 


of  source  dimensions  rather  than  stress  drop  (Thatcher  and  Hanks,  1972).  In 
either  case,  the  lower  the  corner  frequency  at  a  given  Mq  for  earthquakes 
of  a  region,  the  easier  they  are  to  discriminate  from  explosions.  Any  differ¬ 
ences  between  the  LASA  and  NORSAR  data  sets  for  any  earthquake  region  can  be 
attributed  either  to  errors  in  determining  the  long-period  signal  levels 
because  few  of  the  estimates  were  made  with  certainty  or  to  the  assumption 
that  R  is  unity  for  both  stations. 

Attenuation  (t*) 

Relative  mantle  attenuation  was  measured  along  paths  from  each  earthquake 
and  explosion  region  to  LASA  and  NORSAR.  The  measurement  of  attenuation  is 
the  gross  factor  t*,  which  is  travel-time  along  the  ray  path  divided  by  the 
average  quality  factor,  Q. 


t*  was  determined  first  by  normalizing  and  then  summing  the  ground  dis¬ 
placement  (t*  =  0)  spectra  for  each  path.  Figure  20  shows  the  normalized  and 

summed  spectrum  for  paths  from  East  Kazakh  to  LASA.  If  the  assumption  is  made 

_2 

that  the  explosion  source  spectrum  has  an  f  fall  off  at  high  frequencies,  then 


A(f)  ^  S 


irft* 


for  frequencies  above  the  corner  frequency,  where  A(f)  is  the  observed  ground 
displacement  spectrum  (corrected  for  instrument  response),  S  is  the  source 
spectrum  scale  constant,  and  f  is  frequency  in  Hz.  Then  after  taking  the  nat¬ 
ural  log  of  this  equation 


in  [A(f )  ]  +  2  •  £n ( f )  =  -TTt*f  +  £n(S) 

If  In  [A(f)]  +  2*S,n(f)  is  plotted  versus  frequency  for  each  path,  the  slope  of 
the  graph  is  -irt*.  Figure  21  shows  the  result  of  this  process  for  the  East 
Kazakh  to  LASA  path.  The  slope  is  calculated  in  the  frequency  range  1.0  to 
2.5  Hz.  where  the  signal-to-noise  ratio  was  highest,  and  presumably  beyond 
the  corner  frequency,  so  that  the  high-frequency  asymptotic  slope  of  f  ^ 
characterizes  the  source  spectrum.  Plots  such  as  those  in  Figure  17  in  the 
five  Eurasian  reports  showed  that  many  of  the  corner  frequencies  were  above 
1  Hz,  so  that  the  data  may  not  be  well  fitted  by  the  above  equation;  however, 
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FREQUENCY  (Hz) 


Figure  20  Normalized  and  summed  spectra  for  the  paths  from  east.  Kazakh 
I  .ASA  and  to  NORSAK. 
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Figure  21  Ln  A(f)  +  2  ln(f)  frequency  for  the  paths  K.  Kazakh  to  LASA  and 
to  NORSAR. 
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Figure  21  (cont.)  Ln  A(f)  +  2  ln(f)  frequency  for  the  paths  E.  Kazakh 

to  LASA  and  to  NORSAR. 
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if  a  higher  frequency  range  were  used,  then  the  average  spectra  would  be 
seriously  contaminated  with  background  noise.  The  effect  of  using  spectral 
estimates  at  frequencies  below  the  corner  frequency  in  the  fitting  is  to  bias 
the  estimate  of  t*  toward  a  lower  value  than  the  actual  one,  which  would 

explain  why  some  of  our  t*  values  are  unrealistically  close  to  zero  or  even 

_2 

negative.  For  all  events,  source  spectra  with  f  falloff  at  high  frequencies 

were  first  assumed  in  order  to  estimate  t*.  For  the  earthquakes  only,  source 
-3 

spectra  with  f  falloff  at  high  frequencies  were  also  assumed.  The  fitting 
equation  for  this  source  spectral  model  is 

in  [A(f )  ]  +  3 •  in ( f )  -  -irt*f  +  in(S) 

If  the  left-hand  side  of  this  relation  is  plotted  versus  frequency  for  each 

path,  the  slope  of  the  graph  is  again  -irt*;  but  t*  is  now  based  on  an  assumed 

-3  -2 

f  source  spectrum  and  will  be  lower  than  that  obtained  from  the  assumed  f 

source  spectrum. 

-2  -3 

The  t*  values  for  f  and  f  source  models  for  all  paths  are  shown  in 

Table  III.  Note  that  for  three  of  the  five  paths  from  earthquake  regions  to 

-3 

NORSAR,  assumption  of  f  source  models  resulted  in  negative  t*  values,  sug- 
-3 

gesting  that  the  f  source  model  is  inappropriate  for  these  earthquake  regions 

or  that  our  fitting  interval  of  1.0  to  2.5  Hz  was  improper,  as  discussed  above. 

-3 

Although  there  is  a  predominance  of  theoretical  models  that  imply  a  f  source 

-2 

spectrum  rather  than  f  (Geller,  1976),  certain  combinations  of  source-re- 

_2 

ceiver  geometry  can  produce  higher-frequency  ('vf  slope)  signals  from  earth¬ 
quake  models.  Note  also  that  in  each  case  the  t*  value  for  the  path  to  LASA 
is  higher,  implying  more  attenuation,  than  the  t*  value  for  the  path  to  NORSAR, 
which  is  largely  under  shield  areas.  This  result  is  consistent  with  the  basic 
visual  appearance  of  signals  displayed  in  Figure  2  through  6.  Low  t*  values 
for  the  paths  to  NORSAR  (t*  from  0.  to  .34)  are  in  agreement  with  other  t*  re¬ 
sults  for  paths  through  shield  regions  (Der  and  McElfresh,  1977).  All  paths 


Geller,  R. ,  1976.  Scaling  relations  for  earthquake  source  parameters  and 
magnitude.  Bull.  Seism.  Soc,  Am.,  66,  1501. 

Der,  Z.,  and  T.  McElfresh,  1977.  The  relationship  between  anelastic  attenua¬ 
tion  and  regional  amplitude  anomalies  of  short-period  P-wave  in  North 
America,  Bull,  Seism.  Soc.  Am.,  67,  1303-1317. 
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Table  III 


Note:  For  zero  values,  t*  value  calculated  was  negative 
which  is  unrealistic,  so  t*  set  equal  to  zero. 


to  LASA  traverse  the  anomalous  upper  mantle  of  the  Western  United  States;  Der 
and  McElfresh  (1977)  have  shown  that  paths  through  this  region  have  higher  t* 
values  than  paths  entirely  through  shield  regions.  The  paths  from  Kamchatka 
and  MILROW  to  both  LASA  and  NORSAR  are  affected  by  the  dipping  lithospheric 
plate  in  the  source  region. 

The  significant  differences  in  t*  between  East  and  West  Kazakh  is  unex¬ 
plained.  Perhaps  corner  frequencies  are  high  but  variable  between  the  two 
regions  because  of  differences  in  shot  media,  or  perhaps  different  depths 
result  in  differing  effects  of  pP.  This  problem  merits  further  study. 


DISCRIMINATION  ASPECTS 


The  first  part  of  this  study  has  been  aimed  at  uncovering  seismic  source 
characteristics  of  the  events  and  propagation  path  effects  on  the  data  set. 

The  second  part  is  aimed  at  extracting  discriminant  parameters  from  the  data 
and  examining  the  utility  of  these  parameters  in  relation  to  the  characteris¬ 
tics  and  effects  uncovered  in  the  first  part.  In  addition  to  comparing  our 
entire  earthquake  signal  population  to  all  explosion  signals,  earthquakes  and 
explosions  lying  in  similar  tectonic  settings  will  also  be  carefully  compared. 

MsA 

In  the  five  Eurasian  reports,  average  Mg  and  m^  values  were  independently 
determined  for  the  earthquakes  (Table  I)  using  selected  WWSSN  and  HGLP  calcu¬ 
lation);  and  the  ALPA,  LASA,  and  NORSAR  arrays.  Average  Mg  values  were  also 
independently  determined  for  the  explosions  in  Table  II.  The  m^  values  for 
the  explosions  were  taken  from  the  NEIS  list  because  independent  calculation 
of  for  the  explosions  fell  close  to  the  NEIS  values. 

The  following  forms  have  been  used  for  magnitude  computations: 
mb  =  log  (^)  +  B(A) 

Mg  =  log  (|)  +  1.66  log  A  +  0.3 

where 

A  =  one  half  the  peak-to-peak  maximum  record  amplitude  reduced  to  nm 
ground  displacement; 

T  *  period  in  seconds  (restricted  to  17  -  23  sec  for  Mg  calculation); 

A  =  epicentral  distance  in  degrees; 

B  =  Gutenberg-Richter  correction  term  for  surface-focus  P  waves. 

For  the  distances  at  which  amplitudes  are  measured  in  this  report,  there 
are  negligble  differences  between  Mg  values  calculated  using  the  traditional 
Mg  formula  above  and  more  accurate  formulas  with  slopes  near  1.0  derived  by 
several  authors. 


i 


Ringdal's  (1976)  method  of  magnitude  averaging,  which  includes  noise 
measurements  as  upper  bounds,  has  been  used.  With  this  method,  magnitudes  at 
the  individual  stations  are  assumed  to  follow  a  gaussian  distribution.  Within 
this  distribution  some  magnitudes  will  fall  below  station  thresholds  based  on 
the  noise  level,  and  Ringdal's  method  substitutes  a  noise  measurement  at  those 
stations  that  do  not  detect  and  then  computes  the  maximum  likelihood  estimate 
of  magnitude  based  upon  measured  signals  and  noise  estimates.  The  effect  of 
this  procedure  is  to  more  accurately  define  the  magnitude  of  events  not  widely 
recorded;  more  specifically,  the  usual  effect  is  to  lower  the  network  magni¬ 
tude.  Because  for  most  events  in  this  study  few  noise  measurements  were  used, 
the  averages  should  be  very  close  to  results  if  only  measured  signal  amplitudes 
were  used.  As  predicted,  magnitude  averages  for  events  with  many  readings  at 
the  noise  level  were  found  to  be  lower  than  what  would  result  if  only  measured 
signal  amplitudes  were  used. 

The  Mg  -  ro^  plot  is  shown  in  Figure  22.  Where  only  noise  measurements 

were  used  in  the  M  determination,  an  arrow  in  the  figure  represents  the  maxi- 
s 

mum  "Ms"  value.  The  line  Mg  =  m^  -  1.5,  that  Blandford  and  Clark  (1975)  and 
earlier  workers  suggested  as  a  suitable  Eurasian  discrimination  line,  sepa¬ 
rated  the  earthquakes  from  the  explosion  population  up  to  m^  ■v  6.  This  line 
would  also  separate  the  suite  of  Asian  events  studied  by  Dahlman  et  al.  (1974). 
Note  that  in  the  Kamchatka,  Tien  Shan,  and  Pamir  regions  events  were  chosen 
which,  from  the  NEIS  list  or  from  Mg  measurements  at  the  long-period  arrays, 
showed  low  Mg  for  their  m^.  Independent  determination  of  Mg  and  m^  for  those 
events  revealed  that  they  were  not  anomalously  low’  in  Mg  for  their  m^.  In 
summary,  Mg  =  m^  -  1.5  is  a  discrimination  line  for  our  data  set  except  at 
high  magnitude  (MILROW)  where  it  should  bend  upward  to  account  for  the  fact 


Ringdal,  F. ,  1976.  Maximum-likelihood  estimate  of  event  magnitude,  Bull. 
Seism.  Soc.  Am.,  64,  789. 

Dahlman,  0.,  H.  Israelson,  A.  Austegard,  and  G.  Hornstrom,  1974.  Definition 
and  identification  of  seismic  events  in  the  USSR,  Bull.  Seism.  Soc.  Am., 
64,  607. 
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Figure  22  Ms  versus  for  earthquakes  and  explosions  from  all  regions. 


that  m^  is  measured  beyond  the  corner  frequency  (von  Seggern  and  Blandford, 
1972). 

Note,  however  that  the  population  limits  are  separated  by  only  0.5  magni¬ 
tude  units.  Then  an  explosion  shot  array  could  move  into  the  earhquake  pop¬ 
ulation.  In  particular  a  West  Kazakh  shot  array  might  easily  masquerade  as  a 
Pamir  earthquake. 

Corner  Frequencies  and  Long-Period  Spectral  Levels 

Figure  23  is  a  plot  of  |Uq|  versus  corner  frequency  as  seen  at  LASA  and 
NORSAR  for  all  earthquake  and  explosion  data.  Few  long-period  P  phases  were 
recorded  at  the  long-period  arrays  due  to  the  small  magnitudes  of  the  events, 

A 

so  values  of  |Uq|  were  estimated  from  the  spectral  data  of  the  short-period 
vertical  component.  Hanks  and  Thatcher  (1972)  showed  that,  for  a  given  long- 
period  spectral  level,  Amchitka  Island  explosions  have  a  higher  corner  fre¬ 
quency  than  Aleutian  earthquakes.  The  physical  basis  of  this  discriminant  is 
the  smaller  source  time  dimension  of  the  explosion  for  a  given  long-period 
signal  level.  However,  Figure  23  shows  no  separation  between  earthquakes  and 
explosions.  There  is  a  weak  trend  toward  higher  corner  frequencies  for  lower 

A 

|Uq|  for  the  explosions  and  the  earthquakes,  but  it  appears  to  be  steeper  than 
the  theoretical  relation  of  f  3  derived  from  the  fact  that  source  dimension 

A 

scales  as  the  cube-root  of  seismic  moment.  In  Figures  2k  and  25  the  lu  I 

'o' 

values  are  shown  by  region  for  the  LASA  and  NORSAR  data  respectively;  earth¬ 
quake  signals  do  not  separate  from  explosions  in  the  nearest  geographic 
regions. 

One  possible  reason  that  the  spectral  parameters  provided  poor  discrim¬ 
ination  was  that  this  study  used  only  short-period  P  recordings,  while  Hanks 
and  Thatcher  used  long-period  body-wave  recordings.  Another  possible  expla¬ 
nation  is  that,  because  of  the  low  magnitudes  of  events  studied  here,  a  clear 


von  Seggern,  D. ,  and  R.  Blandford,  1972.  Source  time  functions  and  spectra 
for  underground  nuclear  explosions,  Geophys.  J.,  31 ,  83-97. 
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Lgure  24  Long-period  spectral  level  versus  corner  frequency  for  earthquakes 
and  explosions  from  all  regiona  from  LASA  P  recordings. 


all  regions  from  NORSAR  P 


estimate  of  the  asymptotic  value  of  |ll  |  was  not  possible  in  most  cases.  Fin¬ 
ally,  the  lack  of  discrimination  could  stem  from  real  characteristics  of  the 
sources  examined  in  this  study. 


Long-Period  Body-Wave  Excitation 

The  ratios  of  long-period  P  and  S  displacement  to  LR  displacement  was 
examined  as  a  seismic  discriminant  for  some  of  the  events.  For  the  Pamir, 
Baikal,  and  Caucasus  earthquakes,  too  few  long-period  P  or  S  observations 
existed  to  determine  the  ratios  of  these  quantities  to  long-period  LR  ground 
displacement;  all  long-period  P  and  S  observations  from  events  in  these  regions 
were  close  to  the  noise  level,  which  is  expected  for  their  small  magnitudes 
(n^  <  5.3). 

The  Kamchatka  earthquakes  were  generally  of  larger  magnitudes  (up  to  = 
6.0),  and  there  were  more  observations  of  long-period  S  and  P  waves  than  for 
the  other  earthquake  regions  studied.  In  von  Seggern  and  Sobel  (1976),  the 
ratios  SH/LR,  SV/LR  and  P/LR  were  compared  to  distributions  obtained  from  LRSM 
stations  for  NTS  explosions  by  von  Seggern  (1972).  The  median  log(SH/LR)  is 
-0.4  for  the  Kamchatka  earthquakes  and  the  median  log(SV/LR)  is  roughly  -0.6. 
Though  these  values  are  higher  than  those  of  typical  NTS  explosions  recorded 
at  LRSM  stations,  which  lie  at  roughly  -1.2,  the  two  populations  overlap  some¬ 
what.  The  median  log(P/LR)  lies  at  roughly  -0.4,  which  is  higher  than  typical 
explosion  values  which  von  Seggern  found,  suggesting  that  P  and  SH  amplitudes 
are  roughly  equal  for  earthquakes  in  Kamchatka. 

Few  SH/LR,  SV/LR,  and  P/LR  measurements  could  be  made  for  the  Tien  Shan 
earthquakes.  These  ratios  were  generally  larger  than  those  of  NTS  explosions 
in  von  Seggern  (1972),  but  ratios  for  the  earthquake  and  explosion  populations 
again  overlapped  somewhat. 

In  general,  long-period  body-wave  excitation  is  a  weak  discriminant  for 
the  events  studied.  This  is  so  because  of  the  apparent  lack  of  complete  sep¬ 
aration  and  because  the  long-period  body-waves  were  close  to  or  below  the 
noise  level,  expected  for  these  small  magnitude  events,  thus  making  applica¬ 
tion  of  the  discriminant  difficult.  A  remedy  would  be  to  employ  Ringdal's 
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(1976)  method  of  magnitude  estimation  with  some  additional  technique  for  cases 
where  no  signals  are  observed. 

Short-Period  Body-Wave  Excitation 

The  ratio  of  short-period  P  displacement  to  S  displacement  was  also  exam¬ 
ined  as  a  seismic  discriminant.  In  the  Pamir,  Baikal,  and  Caucasus  regions 
there  were  too  few  short-period  S  observations  (mostly  unreliable  with  small 
signal-to-noise)  to  determine  the  ratio  of  S  to  P  ground  displacement. 

In  the  Kamchatka  region  there  were  more  observations  of  short-period  S 
than  for  the  other  earthquake  regions  studied.  In  von  Seggern  and  Sobel 
(1^76)  the  ratios  SV/P  and  SH/P  were  compared  to  distributions  obtained  from 
LRSM  stations  for  NTS  explosions  by  von  Seggern  (1972).  The  spread  of  S/P 
ratios  for  the  Kamchatka  earthquakes  was  over  two  orders  of  magnitude,  reflect¬ 
ing  the  variations  in  S/P  ratios  predicted  by  radiation  patterns  for  double¬ 
couple  models  of  an  earthquake.  While  the  NTS  explosion  population  overlapped 
the  Kamchatka  sample  somewhat,  all  but  one  of  the  nine  earthquakes  had  at 
least  one  log(S/P)  which  exceeded  0.0,  a  value  not  nearly  reached  by  any  of 
the  similar  explosion  log-ratios  in  von  Seggern  (1972). 

One  or  more  SV/P  ground  displacement  ratios  were  available  for  8  of  the 
10  Tien  Shan  events.  The  median  of  the  Tien  Shan  data  lay  at  least  one  order 
of  magnitude  above  the  NTS  explosions  (Sobel  and  von  Seggern,  1976). 

In  general,  short-period  body-wave  excitation  is  a  weak  discriminant  for 
events  studied  here  for  the  same  reasons  that  long-period  body  waves  are  weak 
discriminants.  Again,  lower  thresholds  or  more  sophisticated  magnitude  estima¬ 
tion  techniques  would  render  this  discriminant  more  useful. 

Recent  results  for  regional  distances  on  the  compressional  to  shear  ratio 
(Pn,Pg)/Lg  are  encouraging,  but  no  attempt  was  made  to  reform  these  analyses 
for  these  events. 

Complexities 

The  "complexity"  parameter  (Lambert  et  al.,  1969)  has  been  computed  as 


Lambert,  D.,  D.  von  Seggern,  S.  Alexander,  and  G.  Galat,  1969.  The  LONGSHOT 
experiment,  Vol.  Ill,  Comprehensive  Analysis,  SDL  Report  No.  234,  Tele¬ 
dyne  Geotech,  Alexandria,  Virginia  22314. 
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seen  at  LASA  and  NORSAR  for  all  the  earthquake  and  explosion  data.  Figure  26 
shows  the  complexity  values  versus  m^.  A  high  value  indicated  relatively  high 
coda  energy  (5  to  35  sec  after  P  arrival)  relative  to  the  energy  in  the  first 
5  seconds  of  the  P  wave.  Complexity  numbers  for  the  earthquakes  were  often 
high  due  to  clear  pP  signals.  Complexity  numbers  for  the  explosions  were  low, 
as  would  be  expected,  but  overlapped  the  earthquake  population.  Low  complexity 
numbers  for  some  of  the  earthquakes  were  a  result  either  of  noisy  records  (the 
calculations  include  a  noise  correction)  or  lack  of  a  clear  pP  signal.  In  gen¬ 
eral,  overlapping  complexity  values  for  earthquakes  and  explosions  indicated 
the  unreliability  of  complexity  measurements  to  distinguish  event  type. 

In  Figures  27  and  28  complexity  values  are  shown  by  region  for  the  LASA 
and  NORSAR  data,  respectively.  Complexity  values  for  Kamchatka  earthquakes 
recorded  at  NORSAR  are  low  because  the  pP  phase  for  the  path  to  NORSAR  orgi- 
nates  near  a  nodal  plane  of  the  focal  mechanism.  The  Kamchatka  signals  record¬ 
ed  at  LASA  do  not  separate  from  the  explosion  population  either.  The  closest 
explosion,  MILROW,  has  a  relatively  high  complexity  value,  possibly  due  to  a 
plate  refraction  effect  because  LASA  is  within  the  shadow  zone  defined  by  the 
ray-tracing  study  of  Davies  and  Julian  (1972). 

The  Tien  Shan  earthquakes  do  not  separate  from  the  explosion  population. 

In  fact,  many  of  the  Tien  Shan  vents  recorded  at  LASA  have  complexities  that 
are  among  the  lowest  observed  Li  this  study  either  at  LASA  or  NORSAR.  Clear 
pP  signals  were  not  present  in  most  cases.  The  nearby  Lop  Nor  explosion  falls 
within  the  complexity  range  of  the  Tien  Shan  earthquakes. 

Complexity  numbers  recorded  at  LASA  for  the  Pamir  earthquakes  separate 
from  those  recorded  at  LASA  for  the  nearby  Fast  and  West  Kazakh  explosions. 

This  separation  partially  results  from  the  apparent  pP  signal  on  LASA  beams 
for  most  Pamir  events.  The  complexity  numbers  recorded  at  NORSAR  for  the 
Pamir  earthquakes  are  generally  higher,  but  they  overlap  the  numbers  for  the 
East  and  West  Kazakh  explosions  recorded  at  NORSAR.  This  overlap  can  be 
attributed  to  the  fact  that  NORSAR  beams  for  Pamir  events  show  possible  pP 
signals  for  only  five  of  the  nine  cases. 

Davies,  D, ,  and  B.  Julian,  1972.  A  study  of  short-period  P-wave  signals  from 

LONG-SHOT,  Geophy.  J.,  29,  185. 
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Complexity  versus  mj,  for  earthquakes  and  explosions  from  all 
regions  from  LASA  P  recordings. 


-68- 


COMPLEXITY 


N0R5RP  DRTR 


KAMCHATKA 
TIEN  SHAN 
PAMIRS 
BAIKAL 
CAUCASUS 
ERST  KAZAKH 
WEST  KAZAKH 
LOP  NOP 

BAIKAL  EXPLOSION 


Figure  28  Complexity  versus  for  earthquakes  and  explosions  from  NORSAR 
P  recordings. 


Complexity  numbers  recorded  at  LASA  and  NORSAR  for  the  Baikal  earthquakes 
separate  from  the  nearby  Kazakh  explosions  and  the  Baikal  explosion,  even 
though  only  a  few  pP  phases  were  observed  for  the  Baikal  earthquakes  on  LASA 
and  NORSAR  recordings. 


Complexity  numbers  for  Caucasus  events  recorded  at  LASA  separate  from  the 
complexity  values  for  W.  Kazakh  and  Southwest  USSR  explosions.  However,  com¬ 
plexity  values  for  Caucasus  events  recorded  at  NORSAR  overlap  complexity  values 
for  W.  Kazakh  explosions.  Most  Caucasus  recordings  at  LASA  and  NORSAR  show  ap¬ 
parent  pP  signals.  The  lowest  complexity  value  for  NORSAR  for  Caucasus  events 
was  for  a  possible  subcrustal  event  in  the  Caspian  Sea.  The  low  complexity 
value  for  this  event  likely  stems  from  the  late  arrival  of  pP  and  associated 
crustal  reverberation  signals  beyond  the  30-second  signal  window  used.  Also, 
one  of  the  signals  recorded  at  NORSAR  for  the  West  Kazakh  event  group  has  a 
high  complexity  value  resulting  from  a  clipped  signal.  If  these  two  cases 
are  deleted,  the  complexity  values  for  Caucasus  events  recorded  at  NORSAR 
would  separate  from  the  values  for  Kazakh  explosions  recorded  at  NORSAR. 


Spectral  Ratios 

Spectral  ratios  were  calculated  at  LASA  and  NORSAR  according  to  the 
formula  (Lacoss,  1969) 


SR 


1.95 

1.55 


A(f )df 


//: 


85 

45 


A(f )df 


where  sums  of  the  equivalent  terms  of  the  discrete  Fourier  transform  replaced 
the  integrals  of  the  displacement  spectrum  A(f)  over  frequency.  These  ratios 
are  plotted  versus  in  Figure  29  for  all  the  data  with  no  attenuation  cor¬ 
rection  (t*  =  0).  Spectral  ratios  did  not  appear  to  depend  on  m^,  Mg,  or 
depth.  In  general,  explosions  show  larger  spectral  ratios  (and  more  high 
frequency  energy)  than  earthquakes,  but  the  two  populations  do  not  separate. 
The  spectral  ratios  were  plotted  versus  m^  for  LASA  and  NORSAR  data,  respect¬ 
ively,  in  Figures  30  and  31.  In  Figures  32  and  33  the  LASA  and  NORSAR  data 

-2 

were  corrected  for  attenuation  by  using  t*  values  determined  for  f  source 

-3 

models  for  earthquakes  and  explosions  and  f  source  models  for  earthquakes. 
Bars  connect  earthquake  values  where  there  are  two  different  t*  estimates 
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30  Short-period  P  spectral  ratio  versus  m]-,  for  earthquakes  and 
explosions  from  all  regions  recorded  at  LASA  for  t*=0. 
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Short-period  P  spectral  ratio  versus  for  earthquakes  and 
explosions  from  all  regions  at  NORSAR  for  t*=0. 
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Figure  32 


Short-period  P  spectral  ratio  versus  for  earthquakes  and 
explosions  from  all  regions  [recorded  at  LASA] .  The  t*  values 
are  for  f-^  source  models  for  earthquakes  and  explosions  and 
f“^  source  models  for  earthquakes. 
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Figure  33  Short-period  P  spectral  ratio  versus  mK  for  earthquakes  and 

explosions  from  all  regions  [record  at  NORSAR].  The  t*  values 
are  for  f-^  source  models  for  earthquakes  and  explosions  and 
f~3  source  models  for  earthquakes. 
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from  the  two  source  models.  In  general,  explosions  show  larger  spectral  ratios 
than  earthquakes,  but  the  earthquake  data  does  not  separate  from  explosion 
data  from  nearby  regions.  Lacoss  (1969)  and  von  Seggern  and  Blandford  (1977) 
also  reported  that  spectral  ratios  overlapped  for  Eurasian  events. 

Radiation  Patterns 

Attempts  to  identify  source  mechanisms  for  the  earthquakes  were  largely 
unsuccessful  in  the  five  Eurasian  reports.  The  efforts  were  unsuccessful 
because  determination  of  the  fault  planes  for  earthquakes  whose  magnitudes 
are  below  =  6.0  is  not  generally  possible  with  the  teleseismic  data,  and 
all  our  earthquakes  have  m^  less  than  6.  First  motions  were  determined  when¬ 
ever  possible  from  collected  data,  which  was  mostly  short-period  recordings. 
Unfortunately,  the  short-period  data  were  not  consistent.  Body-waves  for  all 
but  a  few  of  the  presumed  earthquakes  showed  both  distinct  compressional  and 
dilatational  first  motions,  suggesting  that  these  events  were  earthquakes. 
However,  the  overlapping  compressional  and  dilatational  first  motions  on  the 
first  motion  diagrams  indicated  that  this  was  not  a  reliable  discriminant. 

For  body  waves  from  earthquakes  with  thrust  fault  mechanisms,  such  as  in  Kam¬ 
chatka,  lack  of  negative  first  motions  at  teleseismic  distance  would  not  aid 
discrimination. 

In  the  five  Eurasian  reports,  attempts  were  also  made  to  find  LR  radia¬ 
tion  patterns  for  the  earthquakes  by  plotting  the  antilogs  of  the  Mg  values 
for  each  event.  Only  a  few  observed  distributions  of  amplitudes  conformed  to 
any  quadripole  radiation  pattern.  The  scatter  of  amplitudes  was,  in  almost 
every  case,  too  great  to  infer  the  source  type  as  either  the  double-couple  of 
an  earthquake  or  the  uniform  dilatation  of  an  explosion.  The  observed  scatter 
of  data  points  resulted  from  the  small  magnitudes  of  the  events  studied  and 
the  large  epicentral  distances  couple  with  propagation  effects  due  to  varying 
velocity  and  Q  structures  along  the  source-station  paths. 


Lacoss,  R. ,  1969.  A  large-population  LASA  discrimination  experiment,  Techni¬ 
cal  Note  1969-24,  Lincoln  Laboratory,  Lexington,  Massachusetts. 

Von  Seggern,  D.  H.,  and  R.  R.  Blandford,  1977.  Observed  variation  in  the 

spectral  ratio  discriminant  from  short-period  P  waves,  SDAC  Report  No. 
TR-76-12,  Teledyne  Geotech,  Alexandria,  Virginia  22314. 
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All  earthquakes  studied  here  had  apparent  pP  phases  recorded  at  LASA, 
NORSAR,  or  several  of  the  WWSSN  stations.  pP  phases  recorded  at  WWSSN  stations 
agreed  with  the  depths  determined  at  LASA  and  NORSAR,  although  for  crustal 
depth  events  there  is  little  pP  raoveout  at  teleseismic  distances.  Thus,  while 
these  events  could  be  tentatively  identified  as  earthquakes  on  the  basis  of  pP 
observations,  some  identifications  would  be  tenuous  because  explosion  recordings 
occasionally  have  high  coda  with  apparent  pP  phases.  On  the  other  hand,  lack 
of  pP  is  not  a  positive  indication  of  an  explosion.  While  no  attempts  to  en¬ 
hance  pP  observations  were  made  for  these  events,  such  procedures  are  now  under 
study;  and  improved  depth  phase  interpretation  should  be  of  great  value  in 
event  identification. 

Higher-Mode  Observations 

Analysis  of  higher-mode  Rayleigh  waves  was  not  reported  in  the  five 
Eurasian  reports.  Forsyth  (1976)  showed  that  higher-mode  Rayleigh-wave 
arrivals  are  present  on  a  large  portion  of  certain  WWSSN  seismograms 
recorded  from  shallow  Asian  earthquakes  over  relatively  short  continental 
paths.  He  also  devised  a  surface-wave  magnitude  scale  based  upon  higher-mode 
amplitude  which  successfully  identifies  several  Asian  earthquakes  whose  fun¬ 
damental-mode  Rayleigh-wave  and  Love-wave  generation  was  anomalously  low.  The 
NORSAR  and  ALPA  long-period  vertical  beams  have  been  examined  for  higher  modes 
from  the  earthquakes  of  this  study,  but  few  cases  of  clearly  recorded  higher 
modes  were  evident,  especially  at  ALPA,  an  observation  attributable  to  lower 
relative  amplitude  response  of  long-period  array  instruments  to  periods  of  5-15 
seconds  compared  to  that  of  the  WWSSN  instruments  and  also  to  some  signal  loss 
in  beamforming  at  3.6  km/sec  on  the  fundamental  Rayleigh  mode.  Thus,  routine 
examination  of  array  recordings  for  higher  modes  should  be  performed  witli 
special  beamforming  and  filtering. 

In  Figure  34  NORSAR  signals,  1024  sec  long,  are  plotted  from  each  of  the 
five  earthquake  regions  and  from  the  E.  Kazakh  and  S.  W.  Russia  explosion 
sites,  along  with  the  output  of  a  multiple-filter  program  for  determining 


Forsyth,  D.  W.,  1976.  Higher-mode  Rayleigh  waves  as  an  aid  to  seismic  dis¬ 
crimination,  Bull.  Seism.  Soc.  Am.,  66 ,  827-842. 
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group-velocity  curves.  These  signals  were  chosen  because  of  relatively  high 
S/N  ratios,  but  the  others  from  each  region  are  similar,  at  least  in  the 
character  of  the  fundamental  mode.  The  symbols  □  ,  o  and  A  are  plotted 
at  the  first,  second,  and  third  highest  excursions  of  the  envelope  over  each 
of  the  narrow-band  filter  outputs,  thus  allowing  easier  identification  of  more 
than  one  mode.  The  velocity  scale  applies  to  the  contoured  plot,  not  the  time 
series.  Higher  modes  are  apparent  for  the  Tien  Shan,  Pamirs,  and  Baikal  earth¬ 
quakes,  and  possible  for  the  Kamchatka  and  Causcasus  earthquakes.  Paths  to 
NORSAR  from  the  first  three  regions  are  the  most  uniform,  which  accounts  for 
better  recordings  of  higher  modes.  Note  that  because  it  arrives  much  later, 
the  higher  mode  identified  from  the  Tien  Shan  earthquake  is  not  identical  to 
the  mode  from  the  Baikal  or  Pamir  earthquakes.  The  higher  mode  from  the  two 
latter  regions  is  the  first  higher  Rayleigh  mode  that  Cramplin  (1966)  identi¬ 
fied  for  paths  from  Central  Asia  to  Uppsala.  The  Kamchatka  signal  is  complex, 
with  a  probable  multipath  phenomenon  in  the  fundamental  mode  because  the  path 
lies  near  the  continental-oceanic  boundary  in  the  Artie.  The  Caucasus  signal 
is  also  complex,  with  the  presence  of  possibly  several  refracted  or  scattered 
wavefronts.  On  the  other  hand,  fundamental-mode  dispersion  curves  for  both 
explosions  are  very  smooth.  Little  or  no  evidence  of  a  higher-mode  arrival  can 
be  seen  for  the  E.  Kazakh  explosion.  The  fairly  strong  10-14  second  period 
arrival  at  3.9  km/sec  for  the  S.  W.  Russia  explosion  is  a  shear-couple  PL  wave, 
and  the  arrival  in  the  period  range  of  9-14  seconds  and  at  3.9  km/sec  velocity 
is  possibly  the  first  higher  Rayleigh  mode.  However,  the  amplitude  of  this 
mode  relative  to  the  fundamental  at  the  same  period  is  nearly  the  same  on  the 
seismogram.  Since  this  is  contrary  to  theoretical  predictions  for  excitation 
from  a  surface-focus  source  (Forsyth,  1976),  the  difference  might  be  due  to  a 
nearly  non-attenuating  waveguide  for  the  higher  mode,  which  has  its  maximum 
particle  motion  deeper  in  the  crust. 

In  summary,  although  special  processing  could  enhance  its  presence,  high¬ 
er-mode  Rayleigh  is  not  immediately  apparent  on  the  majority  of  seismograms  at 
NORSAR  for  earthquakes  in  this  study.  Also,  higher-mode  arrivals  were  found 


Crampin,  S.,  1966.  Higher  modes  of  seismic  surface  waves:  Propagation  in 
Eurasia,  Bull.  Seism.  Soc.  Am.,  56 ,  1227-1240. 
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for  an  explosion  on  NORSAR  recordings,  indicating  that  the  presence  of  such 
arrivals  is  not  sufficient  to  identify  an  Asian  event  as  an  earthquake.  Stud¬ 
ies  of  seismograms  from  earthquakes  and  nearby  explosions,  such  that  paths 
are  nearly  common,  are  required  to  evaluate  fully  this  possible  discriminant. 

Correlation  with  Active  Faults 

Although  the  location  of  an  event  near  a  known  active  fault  does  not 
imply  that  it  is  necessarily  an  earthquake,  the  location  might  provide  signi¬ 
ficant  diagnostic  information  to  be  weighed  along  with  quantitative  discrim¬ 
ination  data.  The  Kamchatka  earthquakes  studied  here  occured  along  the  coast¬ 
line  within  the  seismic  zone  at  the  top  of  the  lithospheric  plate,  which  is 
dipping  to  the  northwest  under  Kamchatka,  and  are,  therefore,  unassociated  with 
surface  fault  traces.  Landsat  Band  7  satellite  photomosaics  were  made  of  the 
Tien  Shan,  Pamir,  Baikal,  and  Caucasus  regions  for  the  Asian  regional  reports. 
Most  of  the  events  studied  in  those  four  regions  correlated  with  presumed 
active  faults  inferred  from  the  satellite  mosaics. 


STEPWISE  CLASSIFICATION  OF  ASIAN  SEISMIC  EVENTS 


In  this  section  .11  the  discrimination  parameters  are  merged  through  a 
multi-dimensional  approach.  The  approach  used  is  the  stepwise  linear  discrim¬ 
inant  analysis  performed  by  the  program  BMD02M  (Jennrich,  1977).  This  program 
was  also  employed  on  a  data  base  from  Southwestern  United  States  events  by  von 
Seggern  and  Rivers  (1977),  and  a  detailed  explanation  of  the  theory  and  use  of 
the  program  for  an  applied  seismic  discrimination  experiment  can  be  found 
there.  Briefly,  the  program  takes  a  given  number  of  variables  measured  on  two 
or  more  designated  groups,  computes  the  mean  of  each  variable  for  each  group 
and  the  variance  of  each  variable  (assuming  equal  variance  for  all  groups), 
and  then  proceeds  stepwise  by  entering  and  deleting  variables  to  form  the  best 
weighted  linear  combination  of  variables  ("linear  discriminant  function")  for 
each  group  such  that  groups  are  maximally  separated.  Preset  "F-to-enter"  and 
"F-to-remove"  levels  of  the  F-statistic  are  intended  to  control  the  final  sub¬ 
set  of  variables  used  in  group  classification.  However,  these  values  (0.1  and 
.005  respectively)  were  purposely  set  low  so  that  all  variables  were  used.  The 
F-statistic  provides  some  measure  of  each  variable's  ability  to  separate  the 
groups,  with  higher  F  values  implying  more  separation  than  lower  ones.  Finally, 
computed  ji  posteriori  probabilities  for  each  event  are  an  indication  of  the 
degree  of  confidence  in  whatever  classification  was  made  for  it. 

Note  that  results  of  this  experiment  are  completely  valid  only  for  events 
in  this  data  base.  Discriminant  functions  applied  to  the  data  base  here  are 
formed  on  a_  priori  knowledge  of  the  true  classification  of  each  event  in  the 
same  data  base,  and  use  of  such  discriminant  functions  for  new  events  would  be 
made  on  the  assumption  that  values  of  discrimination  parameters  for  the  new 
events  possessed  a  mean-variance  structure  identical  to  that  of  the  training 
set.  This  assumption  would  be  tenuous  because  of  the  small  sample  of  events  in 
the  actual  data  base  used  here.  The  purpose  of  this  experiment  is  firstly,  to 
exploit  a  multi-dimensional  data  set  of  discrimination  parameters  to  determine 
if  improved  classification  might  result  compared  to  simply  using  Mg  -m^,  for 

Jennrich,  R.  I.,  1977.  Stepwise  regression,  in  Statistical  Methods  for  Digi¬ 
tal  Computers,  K.  Enslien  et  al.,  eds.,  John  Wiley  &  Sons,  New  York,  NY. 
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instance,  secondly,  to  quantitatively  focus  on  the  discrimination  capability 
of  each  parameter  measured  for  the  Asian  events;  and  thirdly,  to  identify  those 
events  in  the  data  base  which  are  hardest  to  classify  overall. 

The  six  parameters  used  in  the  stepwise  discrimination  are: 

1)  m^  -  conventional  body-wave  magnitude  average  for  an  event; 

2)  Mg  -  surface-wave  magnitude  for  an  event,  computed  by  Ringdal's 

method  if  possible; 

3)  f  -  corner  frequency  estimated  from  the  attenuation-corrected 

spectra; 

4)  Cp  -  complexity  of  the  P  wave  as  determined  over  the  first  35  sec 

of  signal; 

5)  SR  -  spectral  ratio  from  a  spectrum  corrected  for  attenuation 

-2 

to  an  f  source; 

A 

6)  Uq  -  long-period  spectral  level  estimated  from  the  attenuation- 

corrected  spectra. 

Note  that  fc>  SR,  and  Uq  are  all  derived  quantities,  not  primary  measurements, 

due  to  attenuation  corrections  made  to  the  spectra.  The  t*  values  used  are 

listed  in  Table  III.  Again,  note  that  these  attenuation  corrections  are  not 

precise  and  that  they  may  have  led  to  erroneous  source  spectral  estimates. 

Yet  the  original  Asian  reports  show  that  f  and  Uq  are  relatively  unaffected 

by  the  attenuation  correction.  For  the  earthquake  SR,  correction  to  an  f  ^ 

-3 

or  f  source  does  not  make  a  great  difference  in  the  positions  of  these 
values  relative  to  explosion  values  (see  Figures  29  to  33),  and  they  should  be 
adequate  for  the  analysis.  Note  that  for  the  last  four  of  the  six  parameters, 
only  LASA  and  NORSAR  values  were  computed.  These  values  will  be  used  sepa¬ 
rately  and  as  an  average  quantity  in  classification  experiments. 

Classification  of  Events  According  to  Regions 

While  the  main  effort  is  aimed  toward  discrimination  between  natural 
and  artificial  events,  the  classification  of  the  earthquakes  according  to 
source  region  on  the  basis  of  parameters  measured  for  the  purpose  of  earth¬ 
quake-explosion  discrimination  presents  an  interesting  analysis.  Results 
of  this  analysis  will  indicate  the  degree  of  regional  effects  on  some,  but 
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not  all,  parameters.  Because  of  their  limited  number,  such  analysis  was  not 
attempted  for  explosions. 


For  earthquakes  in  the  five  regions,  the  LASA  and  NORSAR  values  for  para¬ 
meters  3  through  6  were  first  averaged  in  this  experiment  while  values  for 
parameters  1  and  2  (m^  and  Mg)  are  from  Table  I;  results  are  listed  in  Table 
IV  at  the  right.  Only  nine  of  the  forty  earthquakes  were  classified  outside  of 
their  true  group  on  the  basis  of  a  posteriori  probabilities.  The  high  F-stat- 
istic  for  m^  is  artifical,  attributable  to  the  fact  that  events  selected  for 
each  region  spanned  somewhat  different  magnitude  ranges.  Thus,  any  differences 
in  spectral  parameters  f ^  and  Uq  are  connected  to  the  magnitude  range  and 
therefore  may  not  represent  region-dependent  source  differences.  However, 
should  be  relatively  independent  of  magnitude  and  therefore  represents  a  real 
source  effect.  Spectral  ratio  is  seen  to  have  the  least  classification  poten¬ 
tial,  and  this  is  not  unexpected  when  consideration  is  given  to  the  wide  range 
of  this  value  as  measured  at  LASA  and  NORSAR  (Figures  32  and  33) . 

The  experiment  was  repeated  using  LASA  and  NORSAR  estimates  separately 
for  parameters  3  through  6,  and  the  results  are  shown  in  Table  IV.  The  pro¬ 
portion  of  misclassifications  dropped  somewhat  compared  to  results  with  aver¬ 
aged  parameters,  suggesting  that  station  effects  are  an  important  aid  in  earth¬ 
quake  classification. 

The  fact  that  the  earthquakes  separated  fairly  well  into  regions  on  the 
basis  of  the  six  measured  discrimination  parameters  was  a  significant  result 
pointing  to  the  need  for  using  regionalized  discriminant  functions.  The  sign¬ 
ificance  is  exaggerated  in  our  experiment  because  of  the  lack  of  a  planned 
sampling  procedure  for  the  earthquakes;  for  example,  all  the  Kamchatka  earth¬ 
quakes  have  larger  m^  than  the  Pamir  ones. 

Discrimination  of  Earthquakes  and  Explosions  on  a  Regional  Basis 

Past  studies  (e.g.,  Liebermann  and  Pomeroy,  1969,  or  Weichert  and  Basham, 


Liebermann,  R.  C.,  and  P.  W.  Pomeroy,  1969.  Relative  excitation  of  surface 
waves  by  earthquakes  and  underground  explosions,  J.  Geophys.  Res.,  74, 
1575-1590. 

Weichert,  D.  H. ,  and  P.  W.  Basham,  1973.  Deterrence  and  false  alarms  in 
seismic  discrimination,  Bull.  Seism.  Soc.  Am.  63,  1119-1132. 
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Results  for  Classification  of 


89' 


Baikal 


1973)  as  well  as  this  analysis,  suggested  that  discrimination  on  a  regional 
basis  will  be  more  successful  than  when  all  events  are  treated  together. 

Through  this  approach  there  is  no  need  to  determine  regional  correction  fact¬ 
ors,  which  are  often  significant  but  difficult  to  accurately  evaluate.  Follow¬ 
ing  this  reasoning,  earthquakes  in  the  five  regions  will  be.  grouped  with  explo¬ 
sions  in  the  data  base  which  are  in  closest  proximity  to  the  region.  Unfor¬ 
tunately,  the  degree  of  proximity  in  most  of  the  five  cases  is  inadequate  to 
confidently  claim  that  regional  effects  are  truly  common  to  the  explosions 
and  earthquakes.  This  is  especially  true  since  the  explosions,  with  one  or 
two  exceptions,  are  in  aseismic  areas.  The  pairings  were  as  follows  (numbers 
in  parentheses  identify  regions): 

1)  Kamchatka/Milrow  (1/6) 

2)  Tien  Shan/Lop  Nor,  E.  Kazakh  (2/10,7) 

3)  Pamirs/E.  Kazakh,  W.  Kazakh  (3/7,8) 

4)  Baikal/Baikal  explosion  (4/11) 

5)  Caucasus/W.  Kazakh,  S.W.  Russia  (5/8,9) 

The  program  was  executed  three  times  for  each  of  these  pairings:  first  using 
LASA  data  only  for  parameters  3  through  6,  then  NORSAR  data  only  for  them, 
and  finally  using  the  LASA/NORSAR  average  value  (the  LASA  or  NORSAR  value 
alone  was  used  if  both  did  not  have  a  value).  Values  of  parameters  1  and  2 
(m^  and  Mg)  are  from  Tables  I  and  II.  Pertinent  results  are  listed  in  Tables 
V  through  IX  which  contain  two-way  classification  tables  and  F-statistics  for 
the  six  parameters.  Misclassif ications  would  be  identified  by  a  posteriori 
probabilities  <0.5  of  belonging  to  the  correct  source  category,  but  in  fact 
none  of  the  events  considered  were  misclassif ied ,  with  a  probability  of  .975 
being  the  lowest  found  for  any  of  the  events. 

Discrimination  of  Earthquakes  and  Explosions  Without  Regionalization 

Although  regionalized  discrimination  is  preferred,  there  is  some  interest 
in  the  discrimination  of  all  events  in  the  data  base  together.  These  results 
are  listed  in  Table  X.  Using  LASA  or  NORSAR  alone  gave  perfect  separation  of 
the  groups,  although  several  events  had  a  posteriori  probabilities  of  being  in 
the  correct  group  which  were  much  lower  than  in  regionalized  experiments,  as 
low  as  .74.  Using  LASA-NORSAR  averages  for  these  four  parameters  also  resulted 
in  no  misclassif ications.  However,  the  a  posteriori  probability  of  the  741021 
Kamchatka  earthquake  was  a  low  .54.  The  Kamchatka  earthquake  has  low  and 
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very  low  complexity  for  an  earthquake,  but  our  examination  (von  Seggern 
and  Sobel,  1976)  of  more  types  of  data  affirms  that  this  event  is  not  an 
explosion.  While  near  misclassif ication  of  this  event  using  the  combined 
array  data  suggests  that  averaging  discrimination  parameters  over  stations 
may  not  be  worthwhile,  the  use  of  merely  two  stations  as  opposed  to  either 
alone  is  not  sufficient  to  state  this  with  confidence.  Furthermore,  von  Seg¬ 
gern  and  Rivers  (1977)  have  shown  elsewhere  that  discriminant  functions  based 
on  a  six-station  network  average  are  superior  to  those  based  on  single-station 
values  in  separating  earthquakes  and  explosions  in  the  Southwestern  United 
States. 

Table  XI  presents  the  results  of  classification  with  short-period  data 
only,  which  comprises  five  of  the  six  variables  used  above.  The  numerous  mis¬ 
classif  ications  indicate  the  importance  of  long-period  data  (M^  in  this  case) 
in  discriminating  Asian  events.  Note  that  spectral  ratio  is  heavily  weighted 
for  NORSAR,  as  Figure  33  would  suggest,  which  may  be  a  reflection  of  the  good 
high-frequency  recording  of  that  array. 

Value  of  Different  Discrimination  Parameters  for  Asian  Events 

The  set  of  parameters  used  thus  far  in  these  stepwise  discrimination  func 
tions  did  not  include  what  are  thought  to  be  among  the  best  discriminants: 
depth  of  focus,  Love-wave  magnitude,  and  shear-wave  magnitude.  Depth  of  focus 
could  not  be  quantitatively  handled  in  the  program,  and  the  viewpoint  was 
taken  that  depth  was  unknown  and,  therefore,  inconsequential  to  discrimination 
No  Love-wave  or  shear-wave  magnitudes  were  entered  because  of  the  paucity  of 
data  gathered  on  these  phases  in  the  original  five  Asian  regional  reports.  So, 
considering  the  six  parameters  available  in  all  cases.  Table  IV  through  XI 
suggest  that  routinely  measured  and  Mg  parameters  are  essential  data. 

Beyond  those  two,  the  classification  power  of  the  remaining  parameters  is 
marginal,  especially  in  the  experiment  treating  all  events  together  (Table  XI) 
Complexity  weighs  heavily  in  a  few  cases,  but  it  is  not  consistent  between 
LASA  and  NORSAR.  The  spectral  parameters  fc,  U^,  and  SR  also  exhibit  this 
inconsistency.  For  all  parameters  but  spectral  ratio,  there  is  no  clearly 
better  performance  at  one  array  than  at  the  other.  By  using  the  spectral 
ratio,  classification  power  at  NORSAR  seems  definitely  to  be  enhanced. 
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Table  V 


Results  for  Classification  of 
Events  from  Kamchatka/Amchitka 


Classification 

Table 


LASA 


NORSAR 


COMBINED 


9  0 

0  1 


F-Statistics 


LASA 


COMBINED 


5.5 

5.5 

0.2 

0.8 

6.2 

0.1 

■k 

* 

A  .0 

2.5 

2.7 

* 

*  indicates  variables  not  entered,  so  F  <  .01 
Note:  no  NORSAR  data  on  Milrow  explosion 


Table  VI 

Results  for  Classification  of  Events 
From  Tien  Shan/Lop  Nor,  E.  Kazakh 


Classification 

Table 


I  .ASA  NORSAR  COMBINED 


F-Statistics 


f 


Uo 

SR 


LASA 

NORSAR 

COMBINED 

24. 3 

22.  A 

25.0 

16.8 

18.8 

21  .2 

0.3 

3.7 

0.5 

0.1 

2.1 

0.6 

0.1 

7.3 

0.7 

2.0 

3 . 6 

0 .  A 

-92- 


Table  VII 


Results  for  Classification  of  Events 
From  Pamirs/E.  Kazakh,  W.  Kazakh 


Classification 

Table 


LASA  NORSAR  COMBINED 

9  0 

0  10 


7  0 

0  8 


10 

0 

0 

10 

F-Statistics 


% 

l 

uo 

SR 


LASA 

NORSAR 

COMB INI 

49.1 

31.9 

34.1 

10.4 

20.7 

23.6 

0.4 

0.1 

0.1 

5.3 

0.1 

0.1 

1.0 

3.9 

5.4 

3.0 

5.4 

6.3 

Table 

VIII 

Results  for  Classification  of  Events 
From  Baikal/Siberia 


Classification 

Table 


LASA  NORSAR 


COMBINED 

5  0 

0  1 


F-Statistics 


NQRSAk 

8.8 

7.4 

* 

* 

12.0 

1.3 


*  indicates  variables  not  entered,  so  F  <  .01 
Note:  no  LASA  data  on  presumed  Siberian  explosion 


COMBINED 

1.5 

7.4 

* 

* 

26.6 

7.8 
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Table  IX 


Results  for  Classification  of  Events 


From  Caucasus /W. 

Kazakh, 

S.W. 

Russia 

Classification 

LASA 

NORSAR 

COMBINED 

Table 

6  0 

6 

0 

7  0 

0  6 

0 

4 

0  6 

F-Statistlcs 


M„ 


Uo 

SR 


LASA 

10.0 

0.3 

16.8 

3.1 

2.6 

0.4 


NORSAR  COMBINED 


12.1 

17.0 

8.6 

13.7 

2.7 

* 

1.3 

0.2 

2.3 

1.4 

3.9 

1.2 

*  Indicates  variable  not  entered,  so  F  <.01 


Table  X 

Results  for  Classification  of  Events 
From  all  Asian  Regions 


LASA  NORSAR  COMBINED 


Classif ication 

33 

0 

Table 

0 

12 

37 

0 

0 

12 

40 

0 

0 

15 

F-Statistics 


LASA 

NORSAR 

COMB IN 

18.0 

44.4 

84.9 

63.6 

32.8 

21.6 

1.6 

2.1 

1.5 

0.6 

* 

0.1 

* 

0.3 

0.4 

0.2 

5.3 

1.3 

*  indicates  variable  not  entered,  so  F  <  .01 
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Table  XI 


Class if ication 
Table 


F-Statistics  J 


*  M8  not  used 


1 


■] 


Results  for  Classification  of  Events 
From  all  Aslan  Regions  Using  All 
Short-Period  Data 


LASA  NORSAR  COMBINED 

32  5 

4  8 


35 

5 

2 

13 

28 

5 

1 

11 

LASA 

NORSAR 

COMBINED 

“b 

18.0 

5.5 

15.0 

< 

* 

* 

* 

CP 

5.2 

0.4 

0.6 

*1 

0.1 

6.6 

0.1 

A 

Uo 

8.7 

1.3 

11.6 

SR 

0.2 

7.5 

4.6 
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SUMMARY 


Earthquakes  from  Kamchatka,  the  Tien  Shan,  the  Pamirs,  the  Baikal  rift 
zone,  and  the  Caucasus  and  explosions  from  the  Aleutians,  E.  Kazakh,  W.  Kazakh, 
Southwest  Russia,  Lop  Nor,  and  Baikal  were  examined  in  the  seismic  discrimina¬ 
tion  context.  One  of  the  major  points  of  the  investigation  was  to  see  whether 
or  not  seismic  signals  from  those  broad  regions  exhibited  any  uniform  char¬ 
acter,  thereby  indicating  uniform  source-path  factors  for  a  given  region. 
However,  observations  of  short-period  P  waveforms  recorded  both  at  LASA  and 
NORSAR  revealed  few  similarities  in  waveforms  from  any  region  to  either  array. 
Observations  of  long-period  Rayleigh  waves  recorded  at  NORSAR  showed  only  some 
similarity  among  waveforms  from  any  region  to  NORSAR,  except  for  the  explosion 
signals . 

Determining  individual  fault-plane  solutions  on  the  basis  of  first  motion 
was  impossible  for  most  of  the  earthquakes  because  of  the  small  magnitudes  of 
most  events.  The  authors  also  were  unable  to  determine  a  radiation  pattern 
from  the  LR  amplitudes,  probably  a  result  of  varying  propagation  effects. 

Corner  frequencies,  long-period  spectral  levels,  and  seismic  moments  were 
estimated  from  the  short-period  LASA  and  NORSAR  spectra.  The  seismic  moments 
indicate  that  the  earthquakes  are  of  intermediate  stress  drop  (tens  of  bars). 

The  authors  measured  t*,  a  measure  of  mantle  attenuation,  along  the  paths 
from  each  earthquake  and  explosion  region  to  LASA  and  NORSAR.  In  general,  t* 
values  for  the  paths  to  LASA  were  higher  (implying  more  attenuation)  than  the 
t*  values  for  the  paths  to  NORSAR,  which  are  largely  through  shield  areas. 

The  following  discriminants  were  most  successful  in  separating  earthquakes 
from  explosions  in  the  data  set: 

-  Surface-wave  versus  body-wave  magnitude  provided  complete  separation 

of  explosions  and  earthquakes;  the  event  magnitudes  were  based  on 

Ringdal's  (1976)  maximum-likelihood  method  in  most  cases.  One  could 

accept  Mfi  ■  m^  -  1.5  as  a  discrimination  line  for  the  data  set.  In 

the  Kamchatka,  Tien  Shan,  and  Pamir  regions,  events  were  chosen  which, 

from  the  NEIS  list  or  from  measurements  at  NORSAR,  were  of  low  M  for 

s 

their  m^.  Independent  determination  in  this  study  of  Mg  and  for  those 
events  revealed  that  they  were  not  anomalously  low  in  Mg  for  their  m^. 
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Although  the  events  studied  here  separated  on  the  Mg  -  m^  plot  into 
earthquakes  and  explosions,  note  that  there  are  other  areas  in  Asia,  notably 
near  30°N,  95  °E,  where  earthquakes  do  not  clearly  separate  from  explosions 
on  an  Mg  -  plot.  Also,  shot  arrays  could  be  used  to  raise  the  Mg  of  an 
explosion  by  about  0.3  Mg  units  or  more  without  affecting  m^,  placing  it 
close  to  the  earthquake  population  studied  here. 

-  All  earthquakes  in  this  study  had  apparent  pP  signals  recorded  at 
LASA,  NORSAR,  or  several  of  the  WWSSN  stations.  Although  little 
moveout  can  be  detected  at  teleseismic  distances  for  shallow  earth¬ 
quakes,  determination  of  this  phase  is  probably  the  second  best 
available  discriminant.  pP  was  not  identified  for  the  explosions 
but  does  not  function  as  a  positive  discriminant. 

Among  the  unreliable  discriminants  in  this  data  set  are:  corner  frequencies, 
long-period  spectral  levels,  long-period  and  short-period  body-wave  excita¬ 
tion,  complexities,  spectral  ratios,  radiation  patterns,  correlation  with 
active  faults,  and  observations  of  higher  modes.  Still,  each  worked  well  for 
particular  station-region  pairs  and  deserves  consideration  for  regional  dis¬ 
crimination.  It  is  suggested  that  regionalized  discriminants  will  have 
enhanced  value  if  the  size  of  the  regions  is  substantially  smaller  than  those 
of  the  five  studied  here. 

Discrimination  parameters  m^,  Mg,  corner  frequency  (f^),  complexity  (Cp), 
spectral  ratio  (SR) ,  and  long-period  spectral  level  (Uq)  were  merged  in  a 
stepwise  linear  discrimination  analysis  program.  When  parameters  were  averaged 
for  each  earthquake  source  region,  the  earthquakes  separated  fairly  well  into 
regions  on  the  basis  of  the  six  discrimination  parameters,  suggesting  the  need 
for  regionalized  discriminant  functions.  When  events  in  each  of  the  five  earth¬ 
quake  regions  and  the  explosions  in  closest  proximity  to  each  region  were  treated 
separately,  the  earthquakes  separated  from  the  explosions.  For  the  entire  data 
set,  using  all  six  parameters  with  LASA  and  NORSAR  average  values  for  f  ,  C  , 

SR  and  UQ,  all  the  earthquakes  separated  well  from  the  explosions  except  the 
741021  Kamchatka  earthquake,  which  has  low  Mg  and  low  complexity,  and  the 
Baikal  explosion,  which  has  very  little  data  and  no  LASA  record.  Using  LASA 
or  NORSAR  data  alone  gave  perfect  separation  of  the  groups.  Of  course,  if  Mg 
were  enhanced  by  a  shot  array  we  would  expect  that  the  stepwise  linear  discrim¬ 
inant  program  would  still  misclassify  the  event.  Using  only  the  short-period 
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discriminants  in  the  analysis  (that  is,  not  using  M  )  resulted  in  numerous 
misclassifications,  emphasizing  the  importance  of  long-period  data  in  dis¬ 
crimination  of  Asian  events.  The  discrimination  would  be  even  poorer  if  com¬ 
plexity  were  enhanced  in  a  shot  array.  Considering  the  six  discriminants 
analyzed  here,  the  parameters  and  are  essential  data,  but  the  classi¬ 
fication  power  of  the  remaining  four  is  marginal. 

The  low  magnitudes  of  the  events  in  this  study  (m^  from  4.2  and  Mg  from 
3.2)  made  clear  discrimination  difficult  with  the  available  data  and  probably 
represent  nearly  the  lowest  threshold  of  multi-discriminant  analysis  for  the 
stations  studied  here.  The  installation  of  high-quality  SRO  statons  in  Asia 
will  permit  the  examination  of  lower-magnitude  events  with  multi-station  digi¬ 
tal  data. 
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